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We present a mesoscopic description of the anomalous transport and reactions of particles in spiny 
dendrites. As a starting point we use two-state Markovian model with the transition probabihties 
depending on residence time variable. The main assumption is that the longer a particle survives 
inside spine, the smaller becomes the transition probability from spine to dendrite. We extend a 
linear model presented in [PRL, 101, 218102 (2008)] and derive the nonlinear Master equations for 
the average densities of particles inside spines and parent dendrite by eliminating residence time 
variable. We show that the flux of particles between spines and parent dendrite is not local in time 
and space. In particular the average flux of particles from a population of spines through spines 
necks into parent dendrite depends on chemical reactions in spines. This memory effect means that 
one can not separate the exchange flux of particles and the chemical reactions inside spines. This 
phenomenon does not exist in the Markovian case. The flux of particles from dendrite to spines 
is found to depend on the transport process inside dendrite. We show that if the particles inside 
a dendrite have constant velocity, the mean particle's position {x{t)) increases as t^ with fj, < 1 
(anomalous convection). We derive a fractional convection-diffusion equation for the total density 
of particles. 

PACS numbers: 



I. INTRODUCTION 



Dendritic spines play a very important role in regulating neuronal activity of the cerebellar cortex because the 
majority of excitatory synapses are located on spines P, [1] . The later are tiny bulbous protrusions from dendrites 
consisting of head (~ 1 fim) and thin neck (~ 0.1 fim). One of the main functions of spines is to help transmit electrical 
signals. The problem of propagation of action potential in the spiny dendrites has attracted enormous attention in 
past years. Several cable models have been suggested in the literature to study the spine-dendrite interaction on the 
macroscopic level. Baer and Rinzel Q developed a phenomenological cable theory for spiny dendrites and found that 
the propagation rate of local excitation strongly depends on a spine-stem resistance. The dynamic structure of spines 
(the changes in the shape and size) has been studied in Q . Modified FitzHugh-Nagumo model jH was the subject of 
research by Coombes and Bressloff. Population of spines can be treated as a continuous function of spatial position 
along the dendrite [1, 0| • It should be noted that these cable models are phenomenological and not derived from the 
electro-diffusion equations for ions in spiny dendrites. 

Recent experiments together with numerical simulations by Santamaria et al @ showed that the transport of inert 
particles in spiny dendrites is anomalous. It was found that the mean-square displacement (a;^(f)) is proportional to 

with /i < 1. It turns out the dendritic spines act as traps of particles. The narrow spine neck significantly decreases 
the effective diffusion of particles into and out of dendritic spines. This results in slow anomalous diffusion along the 
shaft of dendrite Q. Based on these experiments, Henry et al [l3| suggested the fractional Nernst-Planck equations 
for electro-diffusion of ions in spiny dendrites and developed a cable model involving time-fractional derivatives. In 
particular, the average fiow of ions between the shaft of dendrite and the spines is not local in time. Note that the 
coupling between spines and dendrites has been studied either phenomenologically or on microscopic level of a single 
spine |lll4l3l | . The influence of a dendritic spine on the spread of calcium in the parent dendrite was studied in [l4| . 
The passive diffusion in a tube with dead ends has been considered in [l^. Recently mesoscopic non-Markovian 
model for spines-dendrite interaction has been developed [igI [l7j . The aim of this paper is t o g ive an alternative 
Markovian model involving the residence time variable. This paper is an essential extension of [16j with new results 
and examples. 

We are concerned with a mesoscopic description of the anomalous transport and reactions of particles in spiny 
dendrites. The mesoscopic approach involves a detailed description of the behavior of particles on the microscopic 
level. At the same time within mesoscopic approach we can introduce the mean densities of particles and neglect 
the random fluctuations around the mean behavior [l3| . The system of linear Master equations for the mean density 
of particles inside a dendrite, ni(a;,t), and the density of particles inside a population of spines, n2{x,t) has been 
derived in [iGj . The starting point was non-Markovian continuous time random walk (CTRW) model involving the 
integral balance equations for average densities. The aim of this paper is to extend this analysis for non-linear case 
and derive a mesoscopic system of equations from Markovian model. One way to deal with non-Markovian process is 
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to introduce supplementary variables that make the process Markovian [l^, [l^ . Here we introduce the residence time 
variable t which is the time interval between arrival of particle in dendrite or spine and the time of leaving them. The 
particle has a zero residence time when it just arrives inside the spines from dendrite or inside dendrite from spines. 
This idea has been used in [Tsl, [20l-[23l . 



II. TWO-STATE REACTION-TRANSPORT MODEL 



We start with the microscopic two-state model of a random particle's movement in a spiny dendrite. When the 
particle is inside a dendrite, it performs a random walk for a random time Ti (dendrite's residence time) before hitting 
the neck of spine. When the particle hits the neck, it is trapped inside the spine for a random time T2 (spine's residence 
time). After spending the random time T2 inside spine, the particle is released to the parent dendrite through spine 
neck. The particle starts to perform a random walk inside the dendrite until it hits the spine's neck again. Clearly 
the essential feature of this process can be described by two-state random process. The particle can be in one of 
the two states: inside dendrite or inside spines. We assume that the residence times or waiting times Ti and T2 are 
non-negative random variables with probability density functions (pdf 's) ipi (r) and "02 (t) respectively. Note that 
the introduction of the residence time pdf 's can be viewed as the probabilistic method to model the variability in the 
shape of dendritic spines and their density on the parent dendrite. The corresponding survival functions 5'i(T) are 
defined as 

ij,{s)ds i^l,2. (1) 
To describe the switching process we introduce two transition probabilities (hazard functions) [Tsj : 

Thus 71 (r)/i represents the conditional probability of transition from dendrite to spine in the small interval (r, t -\- h) 
given that there is no transition up to time r. The product 72(T)ft, has a similar meaning for the transition from spine 

rZ _ 

to dendrite. We also assume that inside spines the irreversible chemical reaction C A- Cb with the rate r2 takes place. 
It describes the removal of the particles by immobile buffers and pumps [ill, [Hi • 

If we assume that the transition probabilities 71 (t) and 72(''') are constants, then we have a classical two-state 
Markovian model. Master equations for the mean density of particles inside a dendrite, ni{x,t), and the density of 
particles inside a population of spines, 712 (x,t), are 

L^m - 7ini -1-72712, (3) 

ot 

_— = (n2)n2 -72^2 4-7ini, (4) 

where the reaction rate r2 (71.2) depends on the local density of particles n2 [l3l- Here is the transport operator 
acting on x-coordinate along the dendrite. 

To define the transport operator L^, one can use the continuous-time random walk (CTRW) model in which the 
particles follow the path of the compound Poisson process. The probability of a jump during time interval of length 
h is X{x)h + o(h) and jumps density function is w(z) [17]. Then 



LxTT-i = X{x)ni + / X{z)ni{z,t)w{x — z)dz. (5) 

Jr 

If we adopt the Nernst-Planck equation for the particles fiow then 

d{v{x,t)ni) 

^^"^ = d-x + "^9^' 

where D is the diffusivity of the particles. The advection velocity v{x,t) = —^id4>/dx, where fj, is the mobility, and 
the electrostatic potential cf). Note that v{x,t) depends on the total concentration of particles (ions) through the 
Poisson equation. 

Our aim now is to extend the Markovian model ([3]) and (|H) to the two-state semi-Markov model [l3| • Applications 
of the semi-Markov processes to chemical kinetics can be found in [2^. Note that stochastic two-state models 
with nonexponential waiting time distributions occur in many areas of natural sciences. We mention the stochastic 
resonance theory [23|, two-state model for anomalous diffusion 25], two-state gating process for ion channels 
propagation of tumor cells f27| , superdifFusion theory and random walk with memory [28j . 
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III. MARKOVIAN MODEL INVOLVING A RESIDENCE TIME VARIABLE 



The purpose of this section is to formulate the Markovian model for the transport and reactions of particles inside 
spiny dendrites. If the transition probabilities 71 (r) and 72 (t) are not constants, it is convenient to introduce the 
densities of particles depending on r [13] ■ Let ^1 (x, r, i) be the density of particles at point x at time t whose residence 
time inside the dendrite lies in the interval (r, r + dr) . The corresponding density of particles inside a population of 
spines is ^2(2;, t, t). Integration of r, t) over residence time variable r gives the mean densities rii(t, x) at point x 
at time t: 

poo 

ni{x,t) = ^i{x,T,t)dT, 1 = 1,2. (7) 
Jo 

This model is similar to well-known age-structured models in which the population density of individuals depends 
explicitly on the age r [3l|. Of course, we should make a clear distinction between the residence time r since the last 
jump and the residence time of a particle from t = 0. The latter is not considered here. 

The crucial question here whether the particles "remember" how long they have been inside spines or dendrite. 
The experiments \§\ suggest the existence of the memory effects which lead to the anomalous transport of particles 
along the dendrite. To model these effects, we assume that the transition probabilities 7i(T) depend on the residence 
time variable t [l^. Thus the probability of a transition from a dendrite to spines during a small time interval of 
length h is 7i(r)/i -|- o(/i), and the backward transition has the probability 72(t)/i -I- o{h). 

Since the average movement of particles inside the dendrite is governed by the operator L^, the balance equations 
for £_i{x,T,t) and ^2(2;, T,i) are 



For a dendrite 



^ + ^ =ix6-7i(^)6- (8) 

For spines 

^ + ^ = ~l2{r)i2-r^{n2)i2. (9) 

The derivation of ([8|) and ([9|) together with the transport operator ([5]) is given in Appendix A. Note that this Markovian 
model can be easily generalized to include various nonlinear terms. 
Initial conditions are 

e»(x,r,0) =n°(a;)/,(r|x) i = l,2 (10) 

where n^{x) is the initial densities of particles and fi{T\x) is the conditional waiting time distribution for the particles 
inside dendrite (i — 1) and spines (i — 2) at time t — {) with the property fi{T\x)dT = 1 [20I [2l|. 

Boundary conditions at r = : 

For a dendrite 

Ux,0,t) = I ^2{r)^2{x,T,t)dT. (11) 

For spines 

Ux,0,t) = / Jl{T)^l{x,T,t)dT (12) 











These boundary conditions have the following meaning. Particles inside dendrite with the residence time t = at 
the point x are created with the rate 72(1") (see pip ). The density of particles just arriving inside dendrite, £,i{x, 0, i), 
can be found by integration of the product "f2{T)£,2{x,T,t) over all residence times. The density of particles just 
arriving inside spines ^2(2:, 0, t) can be found in a similar way (see (dU). 

Now we assume that the convective velocity v{x, t) = v = const and the jump rate X{x) = const. Using the method 
of characteristics and the Fourier transform, we obtain the solutions to ([5]) and ^ (see Appendix B): 

For a dendrite 

^i{x,T,t)^e~fo'>'^-'^''' f ^i{z,0,t^T)p{x-z,T))dz T<t, (13) 
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For spines 

^2ix,T,t)=^2{x,0,t-T)e-^o^2{s)ds-j;_^r-{n,(x.s))ds ^ ^ (^5^ 

Here p{x,t) is the Green function for the transport equation 

I^L.p (17) 

with the initial condition 

pix,0) = S{x). (18) 
In particular, when the advection velocity v = const in ([6]) we have 



(a; - vty 



(19) 



Note that ([8|) or (|9]) has different solutions for t < t and t > t. When t < t, the solution p3|) describes the evolution 
of particles that arrived inside dendrite from spines after t = 0. The formula gives the evolution of density of 
those particles that were present inside dendrite at i = 0. 

The formulas (fT3)) and (fT4)) can be generalized for the case when the advection velocity v and jump rate A depend 
on space and time. We can write 

a(x,T,i) = e-^^^^(^)'^^- /ei(^,0,t-T)p(a;,T|z))dz t < t, (20) 



^i{x,T,t)^e-^^-*'''^^''^'^'' Ci{z,T-t,0)p{x,t\z))dz T>t, (21) 



where the Green function p{x,t\z) can be interpreted as the probability density function for a particle which starts 
at point z in a dendrite and arrives at point x at time t without trapping in spines up to time t. The pdf p(x,t\z) 
obeys the equation dp/dt = LxP with the initial condition p(x, 0|z) = 5{x — z). Of course for nonhomogeneous case 
such as ([5]) or ([6l) the explicit expression for p{x^ t\z) is not available. 

The density ^i{x,T,t) can be also interpreted as the probability density function of finding the particle inside the 
dendrite (z = 1) or spines (i = 2) at the point x at time t such that the residence time lies in the interval (r, r + dr). 
Since 5'i(r) = ipi{T), it follows from ^ that the survival function ^'i(T) is 

■^,{t) = e-^o'''^'^'^' i^l,2. (22) 



One can see that this exponential factor appears in the solutions p3|) - p6)) . That is why p3| - (|T6l) have a very simple 
probabilistic meaning of the law of total probability. Note that the residence time's density '0i(r) can be written in 
terms of the transition rate 7i(T) as follows [3 

Mr) = l^{r)e- ^0 ^^^''^'"^ i = l,2. (23) 

It is natural to assume that the longer a particle survives inside spine, the smaller becomes the transition probability 
from spine to dendrite. In this case the transition rate 72(7") is a monotonically decreasing function of residence time 
r. For example, if 

= ^ (24) 

P + T 



then the survival function 

«'2(t) 
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/3 + T 
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Thus the assumption of the dependence of transition rate 72 (t) on the residence time like (p4)) leads to a power-law 
probability density function of residence time inside spines 



The experimental evidence to support the hy pot hesis of a power- law distribution like psp with /i < 1 is the subdiffusive 
transport of particles in a spiny dendrite llal • In what follows we will use this distribution to derive the fractional 
equations for the densities of particles. 

IV. NON-MARKOVIAN TWO-STATE MODEL 

The aim of this section is to set up a non-Markovian model for the transport and reactions of particles in spiny 
dendrites by eliminating the residence time variable r. Let us denote the densities of particles just arriving in a 
dendrite and spines at point x at time t by ji{x,t) and i2{x,t) respectively. Obviously 

j^{x,t)^^,{x,Q,t) z = l,2. (26) 

Note that in the paper (l6| ji (a;, t) denotes the number of particles arriving at point x inside dendrite at time t through 
a single spine stem (not a population of spines considered here) and j2{x,t) is the number of particles arriving at 
point a; in a single spine at time t. 

To derive the balance equations for ji (a;, t) andj2(a;,t) we substitute P^ - pB|) into the boundary conditions ([TT|) .p ^ 
and use pO|) . ((23|) . It gives the following equations for ji{x,t) and j2{x,t) 

Ji{x,t) = rV2(r)j2(x,t-T)e--^*-'^="("^(^'^»'^^dT 

JQ 

-fn°(a;)e~^o'-2"(»2(:'^,«))rf« / ^2(t)*^1(t - t)/2(T - t|x)dT, (27) 



j2ix,t) = I I 'ipi{T)ji{z,t- t)p{x - z,T)dzdT 

V'i(r)*j;^(r - t)n°^{z)fi{T - t\z)p{x - z, T)dzdT. (28) 
Balance equations for ni{x,t) and n2{x,t) can be found by substitution of ((T^ -([TB | and ([TU)) into (O 
ni{x,t) ^ J J i{''')ii{z,t — t)p{x ~ z,T))dzdT 

.+ f [ ■^^{T)-^-\T-t)n\{z)h{T-t\z)p(x-z,T)dzdT, (29) 



n2(,x,t) = f\2{T)j2{x,t - T)e- ^'-^'--^'''^''^'^^^'dT 
Jo 

/oo 
*2(r)*^^(T-i)/2(T-t|a;)dT. (30) 

Note that similar equations have been formulated (not derived) in [16]. For nonhomogeneous case such as ^ instead 
of p{x — z, r) we should use p{x, t\z). In what follows we use the (5-function for residence time distribution at t = : 

Mt\x)=S{t) ^ = l,2. (31) 

It corresponds to the case when the residence time of all particles at t — equals to zero. Substitution of (I5T]) into 
(l?7l) - (|30l) and rearrangement of the integration variables, t — t' — t, give 

Mx,t) = f Mt~t')j2{x,t')e-^^'"'^^^-<^-''^^^'dt' 
Jo 

+n°{x)Mt)e~-^°''^^'"'^"''"^^'^", (32) 
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j2ix,t) ^ / Mt - t') Jliz, t')pix - Z,t-t')dzdt' 

Jo Js. 

+'01 (i) J ni{z)p(x — z,t)dz. 
The balance equations for the densities ni{x,t) are 

ni{x,t) ^ [ ( '^i{t-t')ji{z,t')p{x- z,t-t'))dzdt' 
Jo Jr 

/ Jii{z)p{x - z,t)dz, 



(33) 



(34) 



n2{x,t) = / ^2{t~t')j2{x,t')e-i>^^''-^''^'^^'''dt' 
Jo 

+n°{x)-^2{t)e-Io^^i"^i^^''))d'^, 



(35) 



To obtain the nonhnear Master equations for ni(x, t) and n2{x, t) we differentiate the densities given by p4p and psp 
with respect to time t 

^ = f 'fi{Q)ji{z,t)p{x~ z,0)dz- I I ^Ji{t-t')ji{z,t')p{x~ z,t~t'))dzdt' 

+ f f ^l(t~t')jliz.t')^^^^^—4^^—^dzdt'-^Pl(t) f n"^(z)p(x - z,t)dz 
Jo Jk ' dt 7r 

+*,(0/n?(z)MfL^rf, 



and 



dt 



(?^2(a;,^)) 



*2(i - t')j2{x,t')e- (»2(:«;,s))<is^^/ _^ „0(^)^^(^)g-/o'r2 {n2{x^s))ds 



Using ([32 |) -([35 l) together with (fT7|) . (fT8| and ^i(O) = 1, we can rewrite the last system of equations in a compact 
form 



— = L^ni +ji{x,t) -J2[x,t), 



(36) 



dn2 , , , , _ , , 
-g^ = 32[x, t) - ji{x, t) - (na) n2, 



(37) 



where the densities ji{x,t) and j2ix,t) describe the flux of particles between a population of spines and a parent 
dendrite. Now we need to express ji{x,t) and j2(x,t) in terms of n2{x,t) and ni{x,t) respectively. We obtain 
(Appendix C) 



Jiix,t)= / K2it~t')n2ix,t')e-^^'^^<'^-'^^''y^'''dt', 
Jo 



(38) 



j2{x,t) 



where Ki{t) is the memory kernel defined by 



Ki{t - t')p{x - z,t~ t')ni {z, t')dzdt', 



i = 1,2. 



(39) 



(40) 
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The nonlinear Master equations (l36l) and (I37p together with interaction terms (|38)) and p9l) is a generalization of 
a linear system of equations obtained in [16]. In contrast to the classical Markovian model ^ and the spines- 
dendrite interaction terms ji{x,t) and j2{x,t) are not-local in time and space. The density ji{x,t) describes the 
average flux of particles from a population of spines through spine necks into parent dendrite. The characteristic 
feature of this flux is that it depends on chemical reactions through the exponential term e~ -^t' {"■2ix, s))ds ^ means 
that one can not separate the flux of particles from spines to a dendrite from the chemical reactions inside the spines. 
Similar effects have been found for the reaction-transport systems in [21I [27l . [29j . This phenomenon does not exist in 
the Markovian case for which the memory kernel is delta-function 

K,it) ^ -f^Sit) (41) 

since Ki{s) = 7^. In this case it follows from (p8|) . (|39|) and (|41|) that the fluxes of particles are local: 

ji{x,t) = j2n2{x,t), j2{x,t) = '-fini{x,t). (42) 

Note that the flux of particles from dendrite to spines (|39|) depends on the transport process inside dendrite. 



V. ANOMALOUS ADVECTION 



Consider the case when there is no reaction inside the spines: (^2) ~ 0. It has been found ^16.] that if the particle 
inside dendrite performs a Brownian motion and the residence time's PDF ip2{T) behaves like 

V'2(r)^(^) , ^i<l (43) 

as r — )• 00, the mean squared displacement, < x'^{t) >, exhibits a subdiffusive behavior t^. The reason for this 
anomalous diffusion is that the population of spines acts as a trap of particles [8] . 

The aim of this section is to show that if the particles move inside a dendrite with constant velocity v, the mean 
particle's position {x{t)) increases as (anomalous advection). For a constant velocity v, the transport operator 
takes the form 

L.n, = -v^. (44) 

First we find the Laplace transform of {x{t)) 

where n(fc, s) = ni{k, s) -\- n2{k, s) is the Fourier-Laplace transform of the total density of particles n — ni + n2'. 



n{k,s)^ / e'^''-''n{x,t)dtdx. (46) 

Applying the Fourier-Laplace transform to ((36)) - (p9l) with (712) = 0, we obtain 

n\{k) (*i(s - ^(fc)) + ^m^2{s)) + nl{k) (^2{s) + ^2(s)*i(s - ^{k))] 

n{k, s) = ^ = ^ ^. (47) 

1 - V'2(s)V'i(s - '^{k)) 

Consider the particular case when the initial density of particles inside spines is zero {n^{k) = 0) and the initial 
density inside a dendrite is delta-function: ni(a;, 0) = 5{x) {n\{k) = 1). Then 

^i{s~ikv) + Ms-ikv)^2{s) 
n[k,s) = . (48) 

1 — '02(s)'i/'i(s — ikv) 

since (p{k) — ikv for the transport operator defined by (jH]). We assume that the residence time's PDFs iPi{t) is 
exponential: 
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with the Laplace transform 



^,(.) = -lL_ (49) 



and = l/(7i + s)- The Laplace transform i]j2{s) corresponding to (1431) can be approximated by 

Ms)-l-{r2sf, ^i<l (50) 

for small s [9]; T2 is a parameter with miits of time. The mean waiting time < r >= T^p2{T)dT is infinite in this 
case. In the limit s ^ 0, we find from (HHl), (HH) and (ISOl) that 



n(fc,s)^— ; — - — — r. (51) 



By using (P5|) we find 

(^^(■s)) — rnr — p 



V 



Thus the average position of particle is 



<-^^)^- r(i + ,),,< ^^' 



which is sublinear. This anomalous advection reflects the memory effect corresponding to the slow movement of 
particles in spiny dendrites due to the power law of the residence time distribution for spines (|43p . 



VI. EVOLUTION EQUATION FOR TOTAL DENSITY IN THE LONG-TIME LIMIT 

In this section we derive a governing equation for the total density of particles n(x,t) — ni(x,t) + 712(0;, i) in the 
limit t — >■ 00. We consider the case when initially all particles are inside dendrite: n(x, 0) = ni(x) and n2{x) = 0. By 
using V'i(s + ikv)/'^i{s + ikv) = 71, we rearrange (l47l) as 

\ 0/ ^ (p{k)n{k,s) 
sn{k, s) - ni{x) = ^— . (52) 

1 + 71*2(5) 

Inverse Fourier-Laplace transform gives 

dn /■* 

— = / G{t- s)L,j;n{x,s)ds, (53) 
9t Jo 

where the memory kernel G defined by its Laplace transform 

G(,s) = (54) 

1 + 71*2(5) 

Note that this kernel is different from the standard one in CTRW models (|^D)| . 

Standard advection- diffusion equation. First let us consider the case when the residence time PDF for spines is 

V'2(r) = 726-^^" (55) 

with ^'2(5) = 1/(72 + 5). We assume that 

Lp{k) ^ ivk - Dk^ . (56) 
In the long-time limit t ^ 00 (s — >■ 0) , we obtain from (1521) the governing equation for the total density 
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where 

V — V, U = U. 

71 +72 71 + 72 

Note that Eq. (f57| is vahd for any residence times PDF's iPi{t) and V'2(t) with finite mean residence times < n > 
and < T2 > respectively. The leading order behavior of the Laplace transforms tpi (s) and tjj2 (s) as s — ^ is 

^,(s) - 1- < T, > s, z = l,2. (58) 

Then the effective velocity t;*and diffusivity D* are 

^ <n> ^ <^> ^ 

< Tl > + < r2 > < Tl > + < T2 > 

Anomalous advection- diffusion equation. Now we consider the case when the residence time PDF "02 (t) behaves 
like The infinite mean residence time leads to the anomalous transport of particles along spiny dendrite. By 

using (|50p we obtain the Laplace transform of the survival probability 

*2(.) - = 4,. (60) 

s s'- 

Substitution of (l60l) into (l52l) gives 

-fiT^ist'nik, s) ~ s^^in°(a;)) = Lp{k)n{k, s), 

where ^p{k) is defined by (|56p . We apply the Fourier-Laplace transform inversion and obtain the fractional advection- 
diffusion equation 

df'n dn d^n 



where 



^^'n 1 d /"* n(x,T) n\{x) 

— -dr 



dtt' r(i - ^) at Jo {t - tY r(i - ^i)tt^ 

is the Caputo fractional derivative [13]. Note that the fractional Nernst-Planck equations for ions' movement have been 
suggested in lOi]. In general, the velocity of particles v is not constant since the electric force acting on particles (ions) 
depends on the density of particles. This dependence is described by the Poisson-Nernst-Planck (PNP) equations. 
The challenge is to derive the fractional analog of these equations from the mesoscopic model. Obviously the fractional 
equation (j6ip is just a first step in this direction. 

In fact our model could be used for other applications involving two states: active phase and a quiescent (immobile) 
phase. The classical Markovian switching models involve the Poisson processes with exponential distribution of 
residence times (see, for example, 32]). The non-Markovian switching model developed here can be useful in 
biological modeling [32|-[35|. front propagation in systems with aging and random switching in velocities (36j . large- 
scale transport of solutes in fractured rock involving mobile /immobile transport with power law memory functions 
1371. directed intermittent search on a tree network [3811. 



VII. CONCLUSIONS 



The aim was to give a mesoscopic description of the anomalous transport and reactions of particles in spiny 
dendrites. We extended two-state linear model presented in 16] and derived the nonlinear Master equations for 
the average densities of particles inside spines and a parent dendrite. As a starting point we used the Markovian 
model with an assumption that the transition probabilities depend on the residence time variable. Motivated by 
the experiments Q on anomalous transport of particles along the dendrite we assumed that the longer a particle 
survives inside a spine, the smaller becomes the transition probability from spine to dendrite. By using power-law 
residence time distributions for spines we found that if particles move inside a dendrite with constant velocity v, the 
mean particle's position {x{t)) increases as with /i < 1 (anomalous subadvection). Fractional advection-diffusion 
equation for the total density of particles was derived. We showed that the interaction terms describing the fiux of 
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particles between spines and parent dendrite are not-local in time and space. In particular the average flux of particles 
from a population of spines through spine necks into parent dendrite depends on chemical reactions in spines. It means 
that one can not separate the flux of particles from spines to dendrite from the chemical reactions. This phenomenon 
does not exist in the Markovian case. The flux of particles from dendrite to spines is found to depend on the transport 
process inside dendrite. These results might have a significant implication for nerve cell signalling. The main reason 
for this is that the effective electro-diffusion of ions inside spiny dendrite can not be separated form reactions inside 
spines. An interesting issue is whether this effect has any significant effect on the transport of Ca^+. One can conclude 

that the exponential factor exp ^— J*, (712 (x, s)) dsj in (p8| explains the effect of limited diffusion of Ca^+ along 

dendrites observed in experiments . The derivation of anomalous cable theory from the nonlinear Master equations 
(p6| and (|37|) together with interaction terms (pS)) and ([39|) will be the subject of future work. 
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VIII. APPENDIX A 



To derive the system of equations ^ and (jH]), we start with the balance of particles at the point x 
For a dendrite 

£,i{x,T + h,t + h) ^ {l~ ^i{T)h- \{x)h)S,i{x,T,t) +h / X{z)^2{z,T,t)w{x - z)dz + o{h). (62) 

For spines 

^2{x,T + h,t + h) = {l-y2{T)h)£,2{x,T,t)-r2 {n2)^2{x,T,t)h + o{h). (63) 

The first equation states that the density of particles {x,t + h,t + h) at point x with the residence time t + h inside 
dendrite at time t + h is the sum of the density of particles with the residence time t at time t multiplied by the 
survival probability 1 — 71 (t)/i — X(x)h and the density of particles that jump from different positions z. The jump 
length x — z is distributed according to the dispersal kernel or jump length PDF w{x — z). The second equation 
describes the balance of particles inside spines. The factor 1 — j2{T)h in the RHS of ([55]) is the probability that the 
particles make no transition from spines to dendrite during small time period (t, t + h]. The last term describes 
the decrease of the density ^2 due to the chemical reaction with the rate r2 (712) that depends of the local density 
of particles 712- The advantage to have a Markovian model is that the balance of particles during time {t,t + h] is 
independent of what happened during the previous time interval (0, i]. 

Subtracting ^1 {x, r, t) from both sides of the balance equation and subtracting ^2(2;, r, t) from ([55]) . dividing by 
h, and letting /i — > 0, we obtain the mesoscopic system of reaction-transport equations ^ and ([9]) with the transport 
operator given by ([5]) . 



IX. APPENDIX B 



First we find the solution to ([5]). We denote the Fourier transform of £,i{x,T,t) by ^i{k,T,t). Applying the Fourier 
transform and convolution theorem to ([8]), we obtain 

^ + ^ = vik)i,-^,ir)L (64) 

where (p{k) is the symbol of the pseudo-differential operator L^. It means that the Fourier transform of L^^^i can be 
written as 

T[LMix,T,t)] = ^ik)ii{k,T,t). (65) 

In particular, for the advection-diffusion operator ([B]) with the constant advection velocity u, we have 

if{k):^ivk~Dk^ (66) 

For the operator ([5]) with the constant transition rate A, we obtain 

if{k)^\{w{k)-l), (67) 

where (p{k) is the characteristic exponent of compound Poisson process with intensity A and w(fc) is the Fourier 
transform of jump PDF 'w{z) (l7| . 

The characteristics of PDE (IM|) are determined by the equation dr/dt = 1. Along the straight lines (characteristics) 
T{t) — t — to and r(t) = t + tq, a partial differential equation (PDE) (|64p is reduced to 

'^'^'f^''^ =mi.-l.irit))i.. (68) 



This equation has the solution 

|i(fc,r(t),t) =6(fc,0,io)e^^'^^*-*°^-^4-'^(-(^))'^^ for T{t)^t-U 

and 

ei(fc,r(i),t) =|i(fc,To,0)e'^W*-^'^^(^(^»''^ for T{t)^t + T^. 
Using the inverse Fourier transform and the facts that 



ii{k,0,t- T)e'^''''^^ = I ^i{z,0,t-T)p{x- z,T)dz for t <t, 



6(fc,T-t,0)e'^('')* = I ^i{z,T -t,0)p{x- z,t)dz for t > t, 

we obtain ^ and (fTil) : 



T > t, 



a(a;,r,0 = e"^--*^^^'^'^' / ^i(z, r - i, 0)p(x - z, /or 
where is the inverse Fourier transform of 

In particular, i) is the Green function for the Kolmogorov-Feller equation 

dp f 

— = —Xp + X / p{z,t)w(x — z)dz 
at Jr 

with the initial condition 0) = S{x) and tp{k) — X{w{k) — 1). 
The equation for the density ^2(2;, t, t) 



56 , 56 



-72 (-r)6 - ("2)6 



5t dr 

can be solved by the method of characteristics in the same way. One can obtain 



6(a:,r(i),t) =6(a;,0,to)e"^'o^^(^(^»'^'^-4'^^("^(^'^»'^^ for T{t)=t-to 



and 



e2(a;,T(i),i) =e2(a;,ro,0)e-^o72(r(.))d.-/„'r-(n.(x,.))d. ^(i) ^^^^^^ 

These solutions can be rewritten in the forms given in (|15p and psp . 



X. APPENDIX C. 

Multiplying ^ and §5^ by g/o '■2"(»2(a:,s))d5 taking the Laplace transform £ {/}, we obtain 
c{jiix,t)efo-2Mx^s))ds^ = [nO(a;) + /:|j2(2:,t)e^*'-^("^("'"»'^*}] M^), 
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Then 

C { ji (X, t)eio '■2" ("2 \=C (712(2;, t)eio '■2" (n2{x.,s))ds \ .7Q^ 

L J L J ^2(5) 

Inverse Laplace transform gives 

Jo 

where K2{t) is the memory kernel defined by 

Ms) 



4'2(s) 



From ([7T|) . we get 

Now let us find the expression for j2{x,t) in terms of rii (a:, t). Taking the Fourier-Laplace transform of ([33)) and 
(IMl) . we obtain 

j2 (fc, s) = (ji (fc, s) + n\{k) )C { Vi (t)p(fc, t)} , 

^i(fc,s) = (ji(fc,s)+n;(fc))£{*i(t)p(fc,i)}. 
Taking into account (|69|) and shift theorem, we find 

C{Mmk,t)}^M-s~^{k)), C{^i{t)p{k,t)}^^^{s~^{k)). (72) 

We have 

32{k,s) = hi{k,s)i^ — ■^-77^ ^fii{k,s)Ki{s-(p{k)). (73) 
*i(s - ip(k)) 

Inverse Fourier-Laplace transform of ([75)1 gives (15^ . 



